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Abstract 

We present a method designed for computing solutions of infinite dimensional non¬ 
linear operators f{x) = 0 with a tridiagonal dominant linear part. We recast the 
operator equation into an equivalent Newton-like equation x = T{x) = x — Af{x), 
where A is an approximate inverse of the derivative Df{x) at an approximate solution 
X. We present rigorous computer-assisted calculations showing that T is a contraction 
near x, thus yielding the existence of a solution. Since Df{x) does not have an asymp¬ 
totically diagonal dominant structure, the computation of A is not straightforward. 
This paper provides ideas for computing A, and proposes a new rigorous method for 
proving existence of solutions of nonlinear operators with tridiagonal dominant linear 
part. 
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1 Introduction 

Tridiagonal operators naturally arise in the theory of orthogonal polynomials, ordinary 
differential equations (ODEs), continued fractions, numerical analysis of partial differential 
equations (PDEs), integrable systems, quantum mechanics and solid state physics. Some 
differential operators can be represented by infinite tridiagonal matrices acting in sequence 
spaces, as it is the case for instance for differentiation in frequency space of the Hermite 
functions. Other examples come from the study of ODEs like the Mathieu equation, the 
spheroidal wave equation, the Whittaker-Hill equation and the Lame equation. 

While many well-developed methods and efficient algorithms already exist in the litera¬ 
ture for solving linear tridiagonal matrix equations and computing their inverses, our own 
method has a different flavour. We aim at developing a computational method in order to 
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prove, in a mathematically rigorous and constructive sense, existence of solutions to infinite 
dimensional nonlinear equations of the form 

f{x) = C{x) + N{x) = 0, (1) 

where £ is a tridiagonal linear operator and is a nonlinear operator. The domain of the 
operator / is the space of algebraically decaying sequences 




de£ 



II^IU 


def 


sup{|a;fc|w^} < oo 
k>0 


( 2 ) 


where 

g def r Ij k - Oj 

\ k\ k> 1. 

The assumptions on the linear and nonlinear parts of Q are that £ : 17® —17®“®^ and 
: 17® —>■ 17®“®*', for some > s^. Intuitively, this means that the linear part dominates 
the nonlinear part. Since 17®' C 17®^ for si > S 2 , one can see that / maps 17® into 17®“®*. 

General nonlinear operator equations of the form f{x) = 0 defined on the Banach space 
17® arise in the study of bounded solutions of finite and infinite dimensional dynamical 
systems. For instance, x = {xk)k>o may be the infinite sequence of Fourier coefficients of a 
periodic solution of an ODE, a periodic solution of a delay differential equation (DDE) or 
an equilibrium solution of a PDE with Dirichlet, periodic or Neumann boundary conditions. 
The unknown x may also be the infinite sequence of Chebyshev coefficients of a solution of 
a boundary value problem (BVP), the Hermite coefficients of a solution of an ODE defined 
on an unbounded domain, or the Taylor coefficients of the solution of a Cauchy problem. 
In the case when the differential equation is smooth, the decay rate of the coefficients of 
X will be algebraic or even exponential [T]. In the present paper, we chose to solve 0 in 
the weighed £°° Banach space 17® which corresponds to solutions. In order to exploit 
the analyticity of the solutions, we could follow the idea of [5] and solve Q in weighed 
Banach spaces. This choice of space is not considered in the present paper. 

Recently, several attempts to solve f{x) = 0 in 17® have been successful. They belong 
to a field now called rigorous numerics. This field aims at constructing algorithms that 
provide approximate solutions to a given problem, together with precise bounds implying 
the existence of an exact solution in the mathematically rigorous sense. Equilibria of PDEs 
miais], periodic solutions of DDEs [S], fixed points of infinite dimensional maps [3 and 
periodic solutions of ODEs BM have been computed using such methods. 

One popular idea in rigorous numerics is to recast the problem f{x) = 0 as a problem 
of fixed point of a Newton-like equation of the form T(x) = x — Af(x), where A is an 
approximate inverse of Df{x), and a; is a numerical approximation obtained by computing 
a finite dimensional projection of /. In 01110 IIS], the nonlinear equations under study 
have asymptotically diagonal or block-diagonal dominant linear part, which helps a lot in 
the computation of approximate inverses. In contrast, the present work considers problems 
with tridiagonal dominant linear part. To the best of our knowledge, this is the first attempt 
to compute rigorously solutions of such problems. While our proposed approach is designed 
for a specific class of operators (see assumptions Q and 0) , we believe that it can be 
seen as a first step toward rigorously solving more complicated nonlinear operators with 
tridiagonal dominant linear part. 

The paper is organized as follows. In Section^ we present a method enabling to compute 
(with the help of the computer) pseudo-inverses of tridiagonal operators of a certain class. In 
Section]^ we recast the problem f{x) = 0 as a fixed point problem T{x) = x — Af{x), where 
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A is a pseudo-inverse, and we present the rigorous computational method to prove existence 
of fixed points of T. In Section we present an application and finally, in Section we 
conclude by presenting some interesting future directions. 


2 Computing pseudo-inverses of tridiagonal operators 

This Section is devoted to the construction of a pseudo-inverse of a linear operator with 
tridiagonal tail (see §)• We begin this Section by specifying the assumptions that we make 
on the growth of the tridiagonal terms. Then we use an LU-decomposition to formally 
obtain a formula for the pseudo-inverse. Finally, we check that the (formally defined) 
pseudo-inverse has good mapping properties (see Proposition 2.3). 

Given three sequences {Xk)k>o, {pk)k>o, i/3k)k>o and x G 12®, we define the tridiagonal 
linear operator (acting on x) C{x) = {Ck{x))k>o of 0by 


— XkXk—1 fkkXk -f Pk^k-\-l7 ^ ^ 1? 


(3) 


and Co{x) = /iqXo + PoXi- Assume that there exist real numbers sl > 0, 0 < Ci < C 2 and 
an integer ko such that 


V A: > 0, 


Afe 


kk 


Pk 

/ 

1 


5 

/ 

^k 


< C 2 and y k> ko, Ci < 


k-k 


(4) 


Assume further the existence of <5 S ( Oj 2 I ko > 0 such that 


V fc > fco. 


Afe 



kk 

5 

f^k 


< S. 


(5) 


Then, under assumptions Q and 0 , £ defined by 0 is a tridiagonal operator which maps 
n® into 11®“®^. Indeed, if a; G 12®, then 


ll'C(a;)|U-s^ = sup{|£fc(a;)|a;fc ®"-} 

fc >0 


< C2 ( sup{|a;fe_i|wfc}-bsup{|xfc|wfc}-b sup{|a:fe+i|a;fc} 

Vfc>i 


k>0 


k>0 


< 00. 


From now on, assume for the sake of simplicity that sjv = 0, that is the nonlinear 
part V of 0 maps 12® into f2®. Since fl® is an algebra under discrete convolutions when 
s > 1 (e.g. see [siiin]), then any N which is a combination of such convolutions maps 12® 
into 12®. Assume that using a finite dimensional projection : K™ —)■ K™ of 0, we 
computed a numerical approximation x such that f^'^\x) « 0. We identify x G K™ and 
X = (x, 0,0,0, 0,...) G 12®. We then try to construct a ball 


Bs:{r) = x + Bo{r) = a; + {a; G 12® : ||x||s < r} = {a; G 12® : \\x - a;||s < r} 

centered at x and containing a unique solution of 0, by showing that a specific Newton-like 
operator T{x) = x — Af{x) is a contraction on B^(r). This requires the construction of an 
approximate inverse A of Df{x) = C{x) + DN{x). In order to do so, the structures of £(a;) 
and DN{x) need to be understood. From 0 and 0, £(5;) is a tridiagonal operator with 
entries growing to infinity at the rate fc®^. Moreover, since DN{x) maps fl® into 12®, it is a 
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bounded linear operator. As mentioned above, the expectation is that the coefficients of x 
decay fast to zero. This implies that a reasonable approximation of Df{x) is given by 

/ \ 

D 0 

At ^5=' , (6) 

l^m m I 

\ 0 l^in+1 l^m+l j 

with D =' Df^'^\x) for m large enough. We wish to find the inverse of At in terms of Z?, 
(/3fc)fc>m-i, {tJ^k)k>m and {Xk)k>m- We assume therefore that 


At a; = 


x = y, 


(7) 


where x and y are the infinite vectors 



f ^0 \ 


7 ?/0 ^ 


Xi 


yi 

X = 


, y = 



[ ' J 


\ ' ) 


The infinite part of Q writes 


/ flm 13m 0 0 ... 

Xm+l Mm+l /?m+l 0 

0 Xm+2 Mm+2 [3m+2 

V./ 


( Xm ^ 
^m+1 


V / 


7 ym XjYi XiiYi—i ^ 


i/m+l 


V 


( 8 ) 


/ 


We introduce the notations of the book of P.G. Ciarlet (see Theorem 4.3-2 on page 142 
in [U): 

^2 — X'm+lj ^3 — ..., bi — ^2 — Cl — C2 — 

and (<5„)„gN defined by the induction formula 

^0 = 1, ^1 = bi, and 6n = bn 5n-i - Qn Cn-i Sn- 2 , for n > 2. 

Note that only the <5„ are really useful. 

Let us define the tridiagonal operator T by 


(9) 



f 

Cl 

0 

0 ... \ 


rjn def 

02 

b2 

C2 

0 ... 



0 

03 

&3 

C3 



1 • 



. ... y 


For any infinite vector x = (xo, ■ ■ ■, 

Xk, ■ 


we introduce the notation 

def / 

XjT — (^05 ■ • • 5 ^m— 

if 

and 

XI 

def / 

— (^7715 • • 

• ) ■ • ■) 


Using the notation ei = (1,0, 0, 0,0, • • • )^, the system (|^ becomes 

Txj — yi Xm Xm-l^l- 
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From Theorem 4.3-2 in m, we compute an LC7-decomposition of the tridiagonal oper¬ 
ator defined in ([^ as T = LjUi, where 


, dof 

Lj = 


I \ f 

5i 


a2¥ 1 0 ... 


0 1 ... 

V.y 


and Ui = 


^2 


V . 


0 ^ C2 ... 

0 0 ... 


( 10 ) 


... y 


Hence, the system ® becomes Ljzj = yi — Am Xm-i^i combined with Uixi = zj, that is 


f 0 0 

s_ 


a2| 10 ... 

0 aal^ 1 ... 


( Zm ^ 
-^m+l 


... J 


^ Vm ^m^m—1 ^ 


ym+1 


combined with 


0 C2 ... 


V . 


0 0 ... 


V 

y 

\ 


\ 



^ ■2^771 \ 


^m+1 

= 

■2^m+l 

y 

1 ) 


1 / 


( 11 ) 


( 12 ) 


Both infinite systems 0 and ( |I2| ) can be explicitly solvec. 
System 0 leads to 

^m Vm Am^m—li 


and for any fc > 1 


-\-k — 2/m+/c H -Lj <2fc—Z+2 •• ^/c+1 r 1/m+fc—Z “r ( tj (3-2 •• ^/c+l r 


/ = ! 


— i I v -*-/ •• f 

Ok 

which we rewrite with infinite matrix/vectors notations as 

— Lj \yi Am Xm — l^l] — Li yj AmXm — l ^ 75 

where 


(13) 


zi 


Zm \ 


( Vm \ 



ym-\-l 

^m+2 


ym+2 


, yi = 


V ' J 


\ y 


( 


def j _1 

, vi = Lj ei 


1 

-“2 I 

«3 02 3^ 
-04 03 02 3^ 


V 


The second system (12) leads to the infinite sum (for any fc > 0) 


4 C.XJ 
_ I -|\Z 

^m-\-k — c ^m-\-k ' / X r ^/c+l (^k-\-l ^m-\-k-\-l-! 

0k-\-l ^ ^fe+Z+1 

which we also rewrite with infinite matrix/vector notations as 


xj = Ui ^zj. 


(14) 
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Coupling (131 and (14), we end up with 

xi = Uj~^zj = Ui~^[Lj~^yj — XmXm-i vi] = Ui~^Li~^yi — X^Xm-iWi, (15) 

where wj = UJ^vj. Denoting [UL the hrst row of the infinite matrix Ui~^Lj~^ 
and {wi)q the first element of w/, we can rewrite the first line of (15) as 


)^Q Vi XmXtn—l ■ (1^) 

We now investigate the finite part of the linear system ([^, which is given by 

f xo \ / 0 

xi 0 yi 


D 


^m—2 

\ X^n-i J 


or, according to (16) 


\ Pm—l Xm / 


ym — 2 

V Um-l j 


D 


Letting 


( Xo \ 


( ^ \ 


2/0 ^ 

Xi 


0 


2/1 


+ Bm-l 


= 


^m—2 


0 


ym—2 

\ ^m—1 j 


^ (C/j Lj yi ('^/)o ) 


\ ym -1 ) 


0 

0 ... 

0 

0 

0 

0 ... 

0 

0 

V 0 

0 ... 

0 

K)o / 


K = D-B^,_,Xr, 


we consider its inverse K We denote the last column of iL ^ by {K its last row 

by and its last (“south-east”) element by Then we obtain 

= K~^yp - Brn-i\{Ur^Li~^)^^yi^ {K~^)c^_^ 

= Q(iWi),„_,|® I (17) 

using the tensor product notation. The last line of this identity reads 

Xm-l = {K~^)r^_^yF — Bm-l Li~^) ^^yj^ (18) 


Coming back to (15) and using (18), we see that 
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XI 


Uj Lj yi XmXm—l'^I 

Ur^Lr^yi 

-X 


J/F — /3m-l I (C^/ ^Lj ^)^^yi'^{K 
Ul~^Li~^yi — Xm Wl^^{K~^)r^_^ T/i?! 

{K- 

)m—l,m—1 


Wl 


(19) 


= -AmQ'U;/| 0 |(if ^)r„_i|^yF 

+ (^Ul~^Li~^ + Pm-lXm {K~^)m-l,m-l ® | ^ ’ 

Putting together 0 and p^ , we end up with 

( K-^ -/3—1 j ® ( {UT^Li-^) \) \ 

(At)-i= r 1 r 1 ^ 


ur^Lr^ + A 


where 


^'= Pm-lXm{K ^)m-l,m-l |w^/| ® ^{^1 ^Lj 


In order to get an approximate (pseudo) inverse of , we would like to get a numerical 
approximation of K~^. However the definition of K involves {wi)q, which cannot be explic¬ 
itly computed. By definition, wi = U]~^LJ^ei, so using again the computations made in 
this Section, we get 

(^/)o = iUF^^i)o 


_ ^0 , '^0 


i <^o 


1=1 


^i+i 


'^7n-\-l 


^So 

1^1 SiSi+i 

Given a computational parameter L, we define 

L-l 

E 

'-'i 

and 


Cl . . . C; 02 ... ai+i- 


_ det <5o <5g 


Cl . . . C; 02 . . . Oi+i, 


(20) 


K = D-(3m-lXr, 


I 

0 

... 0 

0 \ 

0 

0 

... 0 

0 

0 

0 

... 0 

0 

V 0 

0 

... 0 

w / 
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We now can consider a nume 
approximate (pseudo) inverse of as 


/ 


= 


consider A^ a numerically computed inverse 
;udol inverse of A^ as 

ur^Lr^ + A 


An 


Mm —1 


of K and then define the 

\ 


where 


A = ^ ® | ^)ro ^ ' 


Lemma 2.1. Assume that m> ko and <5 < |. Then Uj ^ maps into 
Proof. Let z/ G fl® and X[ = Uj^zj. Using (14) and the formula above, we get 


I., I ^ I.. 

pm+fe| ^ |r I Ui 
|Ofc+l| 


ri+k\ + 'y ^ ^ I |cfe + l| |cfc+i| |Zm+fe+i 

|5k+!-!-i| \Zm+k+l\- 


/ 141 , 

- |4+r| IWiI 

Now remember that for all k > 2, dk = bk 4-i — ak Ck-idk- 2 , so 

|4| ^ _ |afc| |cfc-i| |4-2| 

Il5fc-ll l^fel ~ l^fcl l^fc-ll 
^^|4-i||4-2| 

14-1 1 ■ 


/ 


( 21 ) 


( 22 ) 


We introduce Uk = 


141 

|4-i| l^fcl 


which then satisfies 


{ Ul = 1, 
Uk>l- 


4 

Uk-l ’ 


V fc > 2. 


The study of the inductive sequence defined as above, but with > replaced by =, yields 
that for any fc, 7 < Wfc < 1 , where ') "= \ + — 5'^ is the largest root of a: = 1 — ^ (see 

Figure . 
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(a) 5 = 0.2 (b) 5 = 0.5 

Figure 1: The iterations of Un+i = 1 — jun with ui = 1. 


We can then rewrite p2|) in order to get 


\Xm+k\ < ^ |5,+i| .. \b,+i\ |z„+,+,| 


oo 


< 


I I , cl 1 

Km+fel + / . 0 


|<5fc+i| .. |(5fe+i+i| 

1 l^fc+/l 


14+ 


Uk+1 Uk+l |4+/+l| 




/-o 

oo 


|4+;| 
jJ |4+/+i| 


< 


1=0 

Ik/ 


1 


Ci7 


ij l\h+i+i 


Ek 


km+fc +/1 


km+fc+/1 


1=0 


7 J {k + I + 1)*^ {m + k + ly 


(23) 


Finally, since S < - < 7 , 


\Xm+k\ {m + < 


k/ll, 1 


(m + k) 


S + SI, 


Ci7 1 - ^ (fc + 1)®^ (m + ky 


and xi e 

Lemma 2.2. Assume that m> ko and <5 < 4 Then LJ^ maps fl'* into 17® 


□ 


Proof. Let yj G 17® and z/ = Lj^yi. Using (131 and the formula above (without the last 
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term since we do not consider here Lj ^{yj — Xm , we get 


^ \^k l\ 

^ |ym+A;| ~t" ^ ^ | —Z+2I |^^fc+l| \ym+k—l\ 


1^1 

k 


< \ym+k \ + ^ it 1 '^ l^fc-i+2| •■ \bk+l\ \y7n+k-l\ 

l=l ' '‘I 

/I I , l^fe-i + ll l^fc-i+2| ■■ \bk+l\ I I 

^ Uk-l+1 Uk |Ofc_;+i| 

where we use the sequence Uk introduced in the previous proof. We get 


1=0 


1 ) l&fc-i+il 


\ym+k — l I 1 


(24) 


and 


\z„i+k \ {rn + kf < 


< 



m + k 
m + k — I 


For any k > m, we then have 

2 *+*"C 2 || 2 /|| 


\zm+k\im + k) < 


C, 


/[ 2 J / C\ I 




1=0 


7/ \ fc + 1 — ^ 


S + Sl 


t fr' 


'=[!] + ! 


7/ \fc + 1 — Z 


S + SI, 


< 


< 




2'*+®^(72 II2/11^ [ 2 ^+^^ , fS\^ k‘‘+^^+^ 


\ 1 - - V7 
\ 7 ^ ' 


which is bounded uniformly in k since the last term goes to 0 when k goes to oo, and the 
proof is complete. □ 

Proposition 2.3. Assume that m > ko and S < ^. Then A maps into 

Proof. Consider y = {yF,yi)'^ € fl®. Let x = {xf,xi)'^ = Ay. Then, by definition of the 
operator A in (21), 


^F — ArjrayF Pm—1 [ l {Ajyi')c^_.^ 




' ro 


^myr ^m—l^(Uj Lj y/j’ (^m)c^_i ■ 
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By the previous lemmas, Uj ^Lj ^yj G and in particular {Uj ^Lj = 

(ur^Lr^yi)^ is well defined and so is xp- 
Using (211 again, 

xi = —Xm 0 yF + Ui~^Lj~^yj+Ayj. 

Remember that wj = Uf^Lj~^ei, so that wj G U® for any s. According to the previous 
lemmas and the definition of A (see (HI)), we see that xj G U®+®^ . □ 


3 Computations of fixed points of the operator T 

Our main motivation for computing approximate inverses is to prove existence, in a mathe¬ 
matically rigorous sense, of a fixed point of the Newton-like operator T in a set centered at 
a numerical approximation x. The Newton-like operator has the form 


T{x) = x- Af{x), 


(25) 


where A is the approximate inverse (21) of Df{x) computed using the theory of S ectio n 
Since / maps fl® into fl®“®^ and A maps fl® into U®+®^ (thanks to Proposition |2.3[ ), we 
see that T maps the Banach space U® into itself. Our goal is to obtain explicit bounds 
allowing us to show that a given T is a contraction on the ball Bx{r), which yields the 
existence of a fixed point of T (and thus of a zero of /). The fixed point theorem that we 


use (see Theorem 3.1) requires bounds on T and its derivative. We get formulas for these 
bounds in Sections 11 and |3.3[ and then explain in Section |3.4| how to use the so-called 
radii polynomials in order to find a radius r > 0 such that T{B^{r)) C B^{r), and such that 
T is a contraction on Bx{r). 

Before proceeding further, we endow U® with the operation of discrete convolution. 
More precisely, given x = {xk)k>o,y = {yk)k>o G we extend x,y symmetrically by 
X = {xk)k(ii,y = {yk)k& where X-k = Xk, y-k = Vk, for A: > 1. The discrete convolution of 
X and y is then denoted hy x * y, and defined by the (infinite) sum 


(x * y)k = 


^ki yk2' 


fc + fc2 — ^ 


It is known that for s > 1, (U®,*) is an algebra (e.g. see [IH]), that is, if x,y G U®, then 


X * y G U®. This will be useful when we shall look for a bound such as (27) below. We start 


with a classical theorem, whose proof is standard (e.g. see the proof of Lemma 3.3 in [2) 
and is a direct consequence of the contraction mapping theorem. 

Theorem 3.1. For a given s > 1, consider T: U® —>■ U® with T = {Tk)k>o> Tk S M. 
Assume that there exists a point x G U® and vectors Y = {Yk}k>o OLnd Z = {Zk(ry\k>o, 
with Yk, Zk{r) G K, satisfying (for all k > 0) 


|(T(x) - x)fc| < Yk, 


and 


sup 

biMeBoir) 


[DT{x + bi)b2]i, <Zk{r). 


(26) 


(27) 


If there exists r > 0 such that ||y -|- Z(r)\\s < r, then the operator T is a contraction in 
Bx{r) and there exists a unique x G Bx{r) such that T{x) = x. 
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We shall see how to get the bounds Y (Section [3.2| ) and the bounds Z{r) (Section [3^, 
and we shall provide an efficient way of finding a radius r > 0 such that ||y -\- Z{r)\\s < r 
(Section |3.4[ ). The first step however consists in looking for bounds on A. More precisely, 
we need some estimates in order to control the action of This is the goal of the 

following Subsection. 


3.1 Some preliminary computations 

We introduce the notations 

0=- and \ (28) 

7 7(1 

Lemma 3.2. Let yi = (j/m, 2/m+i, • ■ infinite vector and xj = LJ^yi. Assume 

that m > ko and (5 < |. Then, for all k > 0, 


iw.i< 

\j=0 iMm+il 


C30 I 

I ai-k \ym+j\ 

■ j=k+l ^ 

Proof. We again introduce zj = LJ^yi. Combining (23) from Lemma 2.1 and (24) from 
Lemma [ 2 .2| we get 




.. 00 fc+Z I I 

^ nk-\-2l-j \ym-\-j\ 

\b,+i\ 


^ 1 = 0 j=0 


l^j+ll 


1 

7 

1 

7 

= y 


/c I I 00 00 I I 00 1 

\ym+j\ nk+2l-j , f)k-\-2l-j 

k k .k-k j 

^V^ |ym+j| 0’‘~^ \ym+j\ 

^.^1 I^J+ll 1- 

I Qk-i 


-k 

\ym+j\ ^ Qj-k \ym+j\ 

Ik-rn+jl Ik-rn+jl 


□ 


' j=k+l ' ' J 

In particular, we immediately obtain the two following corollaries 
assumptions of Lemma 3.2) which will be useful in the sequel. 

Corollary 3.3. " ’ ' ' .n A,t,, 


r, we immediately obtain tne two toiiowmg coronaries (always under the 
Lemma 3.2) which will be useful in the sequel. 

Recall (28). Then, for wj = {wm,Wm+i, ■ ■ ■)'^ '= Lj^ei, we have 

\wm+k\ ^Vd’^rk, forallk>0. 

\ym\ 

Corollary 3.4. If y is such that ym+k = 0 for any k > n, then 

yk<n- 2 , \x^+^\<y(y^ 0 '^-^kk+ 


(29) 


(30) 


and 


V fc > n — 1 
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More generally, we will also need in the next two Subsections a uniform bound on 
\xm+k \ ijn + for k large enough. We assume here that m > 2 (which will always be 

the case in practice), and define for any integer M 


//I 71 r \ det „M. M 

X = x(.S, m, M, s,sl) = d^ — 


■M 


s+sl 


I gVM ^ Qd+dJ. 


1 / m + M 

l-O \m + M-VM-lJ 


Proposition 3.5. Suppose that M satisfies 


M > max 


—min '/O — s — sl — ^ — y {mlny6 + s + sl + 1)^ — 4mIn Vd 4 


2 In 


(ln0) 


2 ’ ‘ 


Then for all k < M, 

\xm+k \ [m + ky^'"^ < 
and for all k > M 


s+sl ^ nk-l 


Cl 


\i=o 


m + I 


1-1 


\Xm+k \ (m + ky^""^ < 


Proof. Thanks to Lemma [T2 


I Wl (m+ <:)■+“ 


\k — l 


< 


Cl 

yhy 

Cl 






k-l 


m + k 
m +1 


m + I 


1-9 


s+sl 00 


(32) 

(33) 

(34) 


E^ 

l=k+l 


l-k 


m +1 


S + SL ' 


+ 


1-61 


Then for k > M, we split the remaining sum 


E^' 

1=0 


_i f m + k 
m + I 


s+si, L 2 J 


= E^ 


k-l 


1=0 


m + k 
m +1 


S + SL k-[Vk]-l 

+ E " 


k-l 


m + k 
m +1 


S + Si 


fc fc fm + k 
< 02 - 


2 \ m 


S + SI, 




m + k 


E 

;=fc-[v^] 

S + Si 


_i f m + k 
m +1 


1 — 9 \m + k — Vk — 1 


< 9 2 


= X- 


MMf m + M \ ^^M^s+sl j_ 


m 


m + M 


1 — 9 \m + M — y/M — 1 


S + Si 


The justification of the last inequality is contained in the following three lemmas. □ 


Lemma 3.6. If M satisfies (32), then for all k > M 

k k ( m + k\ M M f m + M\ 


02 - 


2 V m 


<9^ — 
2 
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Proof. For a; > 0, let i.pi{x) = 9^x{m + whose derivative is 

(p[{x) = Vd x{m + x)®+®^ + (m + x)®'*’®^ + (s + sl)x (to + x)®^®^~^j 

= (to + x)®~''®^~^ Vo (to + x)x + (to + x) + (s + s^x^j 


= (to + x)®"'’®^ ^ Vo x^ + ^TOlnv^ + s + sl + 1 ^ X + TO^ . 

For 0 < 0 < 1 , the discriminant of in VOx"^ + {m in VO + s + sl + 1 ^x + to given by 


A = 


^TOln-\/0 + s + Si + ij — ^mhwfO, 


is positive. Since M satisfies (321, ViV) < 0 for any x > M and so (pi{k) < (pi{M) for all 


k> M. 


□ 


Lemma 3.7. If M satisfies (32), then for all k > M, 

^■\/fc^2®+®t < qVM n^ + SL 

2 - 2 ■ 

Proof. Let f 2 {x) = O'^x. Then 

VV) = 0'^ ^^^x + 1^ = ^ (v^ln6» + 2) . 


Hence, for x > 


(ln0) 


2 ’ 


W) ^ 0 and so ip 2 {k) < ‘P 2 {M) for all k > M. 


□ 


Lemma 3.8. If M satisfies (32), then for all k > M, 

s+sl 


1 


TO + fc 


1 — 9 \-ui + k — Vk — 1 


< 


1 


m + M 


1 — 9 \]ji + M — VM — 1 


S + Si 


Proof. Let ipz{x) mPx —Then 

TO + X — ^Jx — 1 


W) = 


. + x- v^-1-(to + x)(i-^)^ x + 2V^-to 


(to + X — Vx — 1) 2Vx (to + X — Vx — 1)^ 

Hence, for x > to, VsV) — 0 and ifsV) < Lp^{M) for all k > M. 


□ 


Finally, we will need to bound the error made by using w instead of {wi)q for the 
definition (21) of H. 


Lemma 3.9. Assume that L > ko and 6 < V Then 


a2L 


|(u;/)o-w| < 


l^ml (1 - 0 '^)' 


( 35 ) 
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we 


Proof. Using ([^ together with the sequence (ui) introduced in the proof of Lemma 2.1 
get 

OO I C |2 


\{wi)^-w\ < 


l^L 


\Si\\Si+i 


|ci|...|cz| \a 2 \...\ai+i\ 


<Mfs- (1 

^ ^ ' '/ii 'fij • ' 


Ui Ul J \U2 UlJ^i 


< ,—, Ve"' 
W\ 

q2L 

l^ml (1 - 0 ^)' 


□ 


3.2 Computation of the Y bounds 

From now on, we shall assume for the sake of clarity that the nonlinearity TV of / in 0 is a 
polynomial of degree two. The generalization to a polynomial nonlinearity of higher degree 
could be obtained thanks to the use of the estimates developed in [S] in order to bound 
terms like 

{x^ * ... * 

where x^,...,x^ G Bo{r). Moreover, as long as one is interested in problems with non- 
linearities built from elementary functions of mathematical physics (powers, exponential, 
trigonometric functions, rational, Bessel, elliptic integrals, etc.), our method is applicable. 
Indeed, since these nonlinearities are themselves solutions of low order linear or polynomial 
ODEs, they can be appended to the original problem of interest in order to obtain polyno¬ 
mial nonlinearities, albeit in a higher number of variables. This standard trick is explained 
in more details in [12], and is used in [18] to prove existence of periodic solutions in the 
planar circular restricted three body problem. 

With this in mind, we are ready to compute the bound Y appearing in Theorem |3.1| In 
everything that follows, j-l, when applied to vectors or matrices (even infinite dimensional), 
must be understood component-wise. 

The main estimate of this subsection, that is the bound on Y, is presented in the following 
Proposition: 


Proposition 3.10. Consider an integer M such that 

—s 


M > max 


ln 6 » 


— TO, TO — 2 


and define Y = {Yk)k>o component-wise by 


(/(x))^| + |/3^_i|77 


, 1=0 


,;l /(^)L+i 

Ih-m+l I 




(36) 


(37) 


^ m+/c — 


1=0 \h'rn+l\ i^k+1 


E 




Ql-k m-l-i 

V0<fc<TO — 3 

1 Aim-1-/1 



r]9 


k I 


IMti 


(38) 
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V . — 
m+/c — 


/m—2 


+ /3^-l ^g/ l/(^)Uz ] ]^gfe|A.n| 


, 1=0 


1^, 


m+Z I 


IMti 


l/OzOI^^ Vm-2<A:<M, 

P' 


i=0 


and 


Y^+k=Y^+M‘^'^, yk>M. 
wi , 1. 


( 39 ) 


(40) 


Then 

Proof. By definition of T, 


^m+fc 

|T(x) - x| < F. 


|r(a;)-x| = \Af{x)\. 

Note that since we suppose that / is at most quadratic, and since x is constructed in such a 
way that Xk = 0 for all k > m, we get the identity {f{x))^_^_j. = 0 for all /c > m — 1. Thanks 
to 

\{Af{x))^\ < if{x))p\ + \f3m-i\\{Ur^LJ^ im)i)o\ 
so that using (|^ with n = m — 1 and fc = 0, we get 


/ m—2 


\{Af{x))p\ < \A^ {f{x))p\ + \Pm-l\ »7 XI 


i\fi^)L+i 


. 1=0 


\p 


m-\-l I 


(^m)c„_j 


which provides the bound (37). 
Using (21) again, 


\{Af{x))j\ < |Am| ^ (^m)T-„,_i /(*)^’ + Pm-1 (^m)m-l,m-l I f{x)l)^ ^ |u;/| + |t/j J 


so using ([^, (30) and (31) (again with n = m — 1), we get 


\{Af{x)) 


=i.\ < 

m-\-k \ — 


/m-2 


(^m)r„_i f{x)F + Pm-1 (^m)m-l,m-l ^ I ^ ^ 


Mix)\ 


771+Z 


, z=o 


\p 


■m+Z I 


rjO 


k l^m| 


I Mil 


„fc_i l/(2;)lm+Z „;_fc l/(2;)lm+; v! n ^ 1 ^ o 

n y 0 —-,-h P > 0 —-—, v0<k<m — 6, 

to it+i ■ ■ 


which provides the bound (|38|), and 


\{Af{x)) 


=i.\ < 

m-\-k \ — 


/m-2 


{Am)r.^_i Y Pm-l (^m)„i-l,m-l M I ^ ' 




m+Z 


, Z=0 


\p 


■m+Z I 


rjO 


k |i^in| 


I Mil 


Vfc>m-2, 


^ 0^ l/i 
z=o 


m+Z I 


which provides the bound (39). Finally, by (36), 6l*(m + ky < 9^{m + M)® for all k > M, 
and we obtain the bound (40). □ 


We present in Section 3.4 the rationale behind the definition of Y^+k for k > M. 


V(+7 
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3.3 Computation of the Z bounds 

In order to compute the Z bounds from Theorem |3.1[ we need to estimate the quantity 
DT {x + y) z = {I — ADf {x + y)) z = {l — z — A {Df {x + y) — z 


for all y,z G BQ{r). We are going to bound each term separately in the next two Sub¬ 
subsections. We introduce the notation 




def 



(41) 


3.3.1 Estimates for (I — AA^)z 

In this Sub-subsection, we present the bound on {I — AA^)z, which constitutes the first part 
of a bound for Z. 


Proposition 3.11. Let M be an integer satisfying ( |36[ ). We define Z^ = {Zl)k>o component¬ 
wise by 

I, I 


ryl 

Z_p — 


yl 

^m+k 


and 


I - A^ K 


I - A^,K 
Wl 


|/3„i||A™|02i 


(42) 


lMm|w^-l(l-0^) 


\Am\m-i,m-i]vO'^^-^r, V 0 < < M, 


\a 


^m+k — ^m+M i J ^ H > M. 
^m+k 


(43) 

(44) 


Then for all z G Bg^r), 

Proof. Thanks to ([^ and (21), 

(AAX)^ = Ajn I Dzf -f 


\{l-AA^)z\ < Z\ 


— fim-l {Uj ^ {Tzi Xm^m-l^l)) Q (^m)c„_i 




and so 


— AjYiDzp -|- /3 

m—l^m —1 ^m — 1 H“ ^m^m—1 

= AjyiKZ p -\- ^'m — 1 (^ 111 ) 0 ^^ 1-1 ’ 

((/ - AX) z)p = {i - AmK^ Zp + [w - (w/)q) Zjn-1 (^m)c„,_i ■ 


For z € Bo{r) we have, using (35), 


/ - A^,K 


^rs , |/3m-l| |Am| 


which provides the bound (42). 
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Using again Q and (21), we get 

/ 


(^AA^ z) J — —Xm 


Dzp 


V 



Wj + (Uj ^Lj ^ + A) {Tzi + XmZm-l^l) 


= Zl + XmWl 

^ /^m—1 —1 Zm — 1 Prn—l (z/ + XjnZm — l^^l)^ 

— Zj A Xm ^ ^Zp A Zm—1 “1“ (Xm—lXm {Am^m—l.m — l Z'm—l (^/)o^ 

= Zj A Xm ^Zm—1 ^^F “t” i^m—lXm {Am)m—l^m—l Z'm—l (^/)o)^ 


((' 


— Zl Xm ( ( I AmKj Zp + Pm—lXm iAm)m — l^m—l Zm—1 (^/ )q) ) 1 


and so 


{(-^ Z^ j — Xm l^m—lXm (^-^m')m—l,m—l ^m-l (jF {'^/)o)^ Wp 


For z S Bo{r) we have, using (29) and (35), 


|(/-AA^)z| 


m-\-k — 


/ - Am K 


\l3m-l\\Xm\0^^ ,, , \^nk\^rn 


Wp + 


\Am\m-l,m-l ) V fc > 0, 


which gives (43), as well as (44) thanks to (36). 


□ 


3.3.2 Estimates for A {Df {x + y) — A'l) z 

This Sub-subsection is devoted to the exposition of a bound for A {^Df {x + y) — A'^) z, 
which constitutes the second (and last) part of a bound for Z. This bound is detailed in 
Proposition |3A8| 

Recall the assumption that the nonlinear part N is polynomial of degree 2. Hence, 
Df (x + y) can be written as a finite Taylor expansion 

Df {x + y)= Df {x) + D'^f (x) (y). 


and 

{Df (x + y)- A'l') z = {Df (x) - A''') z -f D^f (x) (y, z). (45) 

We are going to bound the two terms of ( |4^ separately. Let us denote by a the coefficient 
of degree 2 of /, that is D'^f (x) {y,z) = 2a{y * z). We bound this convolution product 
thanks to the following result: 

Lemma 3.12. Let s > 2 be an algebraie decay rate and n > 6, let L > 1 be computational 
parameters. For x,y G W and for any k > 0, 


\{x*y)yf\ < al{nY 


.bll 
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where 




1 + 2 ^ 7 ^ 




fc = 0, 


2 + 2Ey^ 




k-1 

E 


2 + 2^ 


■ + 2 


1 < /c < n, 


4 ln(n — 2) tt^ — 6 


(s — l)L®“i \n —1^ 

Proof. See [ni for a proof of this bound and m for a similar bound for 1 < s < 2. □ 

Remark 3.13. It is important to notice here that ot’^in) = a-’fin) for all k > n. From now 
on, we assume that m is taken larger or equal to 6, which will allow us to use Lemma 3.12 
with n = m. Note that this condition is not stringent, since in practice more than 6 modes 
are usually needed in order to get a good numerical solution x. 


We begin by bounding the first term of ( |45| ). 

Proposition 3.14. Define = C^{x) = (C'fe(S))j,>o component-wise by 


Clix)^^2\a\ 


\xi\ 


l—m—k 


, Vl<fc<m — 1, 


and 


CUki^) = 


2|cr|a^(m) ||x|| 


V fc > 0. 




m+/c 


Then for all z G BqIj) 

I [Df{x) — A'^) z\ < C^{x)r. 
Proof. According to the definition of in ([^, we see that 


{{Df{x) - At) z)^ = {Df{x)z)p - Df^'^\x)zF - 


— 2a ({x * z)f — {x * zf)f) 


( 


y jdm—l ^ra j 


k> n. 


where in the convolution product, Zf must be understood as the infinite vector {zf, 0 ,..., 0 ,.. 
Therefore, ((Zl/(x) — At) z)^ = 0, and for all z G Bo^r), 

\{Df{x) - A^) z\^<2\a\r ^ V l<fc<m-l. 

l—m—k 

Then, remembering that Df{x) = C DN{x) and , we see that 
{{Df{x) - At) z)j = {DN{x)z)j = 2a{x* z)j , 
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so that using Lemma 3.12 for all z G B^ir), we end up with the bound 


!(!,/(.)v.> 0 . 

^m+k 


□ 


We now bound the second term of (451. 

Proposition 3.15. Recall (41) and define = {pk')k>o component-wise by 


UJt 


Then for all y,z G B^ir) 


\D^f{x) (y,z)| <CV. 


Proof. Remembering that D^f (x) {y,z) = 2a{y * z), this is a consequence of Lemma 3.12 

□ 


Finally, 

IA {Df {x + y)- A'f) z| < |y4| {C\x)r + , 

and we are left to bound |A| C^(x) and |A| C^. 

Proposition 3.16. Let M be an integer satisfying (32) and (36). We define 
component-wise by 


k>0 


D^f 


(X) - |A™| CAx) + 


(46) 


.l det / I . I 2 |/3m-l| ^ I'^l \^Tn\ 


Dl^kiS:)=[\AAr^_^CAx) + 


Ci(l-0)a;^^ 


S + Si 


Ih-T. 


2 ry|cr|a^(m)||a;||s ffk-i f ^ 

\m + l) 


S + SI, 


e 


+ T30 ’ V0<A:<M, (47) 


Dl+Mix)=[\AAr^_^CUx) + 

, 2??|cr|a^(m)||x|| 


2 |^m-l| \Am\m-l,m-10W\(^L(.m)\\x\\s\ \Xm\nM 

' 7]t ^ 


Ciil-0)AJ 


S + Sl 




1-0 r 


and 


Then 


Dl^k (S) = dI+m (S) , V fc > M. 

^m+k 


\A\C^ix) < D\A- 

Proof. Thanks to (HI), 

(|A| C\x))f < l^ml CU^) + l/3m-il \Ur^LJ^C}{x)l |A 


ITti 


(48) 


( 49 ) 


0 I Cm — 1 ’ 
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and using (33) 






, < 
s+si, — 


27?|g|a^(m)||a:||, 


so that Mhl) holds. Still thanks to (21), 


(1^1 C\x))j < |A^| CUx) k/| + \Ur^LJ^C}ix)\ + |A™| |/3„_i| \Ur^LJ^C}ix)l |a;z| 


— l-^ml ( + 


Using (29) and (33), we get 


(|>1| C^{x))^^^ < ( |^m|r„_i Cp{x) + 


2 |/3^i |V \<j\ g^Mllxll 

2 \Pm-l\ |^mL_i,^_i V kl amMPII. 


\wj\ + \UY^L-J^C]{x)\. 




S + Si 


JA. 


ml nk 


I Mr? 


+ 


27?|g|a^(m)||a;|| 


E 

\l=0 


3fc_i W + fc\ 

+ l 1 


s-\-sl 


m ■ 


1-6 


, yo<k<M, 


so that (47) holds, and using (29) and ( |34[ ), we get 

2 |/3™_i| rj\a\ a^^{m)\\x\\ 


{\A\C\x)) 


m+M 


< |A 


Cp{x) + 




,\^rn\ 

'l- 1 ^ 

\}^m\ 


+ 


2 ? 7 |g|a^(m)||x|| 


X + 


0 


1-9 


so that (48) holds. As before, (49) follows from (36). 


□ 


We get similar results for the second order term. 

Proposition 3.17. Let M be an integer satisfying (32) and (36). Define 
component-wise by 


= |A^|C| 


2|/3m-i|?7k|a™(m) , . 


D 


2 

m+fc 



2 l^m-ll \AmL_p^_^ 7? |cr| a^(m) 




S + Si 


1 


I'^ml nk 

iMmI 




(E"' 

\i =0 


m k 
m +1 


S + Sl 


1 - 0 


y o<k<M, 


n 2 


and 


Then 


= ( |A. 




2\p^_i\\A,n\^_p^_^p\a\af^{m)\ \X^\ 2p\a\a^^{m) f 6 


D 


m-\-k 


jj2 


Ciil-0)uj^^^ 

^m+M 


IM' 




X + 


1-1 


m+M s 


y k> M. 


^m+k 

'i2 ^ 7 -i 2 


\A\C^ < DT 
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Finally we can sum up all the computations of this Sub-subsection and state the following 
result: 


Proposition 3.18. Let M be an integer satisfying (32| and (36). We define (resp. D^) 
as in Proposition 3.16 (resp. Proposition 3.17) and let 


Z\r) = D\x)r + DW 


Then for all y,z G B^ir) 


A (Df {x + y) — A^') z < Z‘^{r). 

Putting this together with Proposition |3.11[ we end up with the following result: 


Proposition 3.19. Let M be an integer satisfying (32) and (36). Let 

Z{r) =Z\r) + Z^{r). 


Then for all y,z G Bo(r), 


\DT {x + y) z\ < Z{r). 


3.4 The radii polynomials and interval arithmetics 

All the work done up to now in Sections and can be summarized in the following 
statement: 


Theorem 3.20. Let s > 1, and s^ > 0. Assume that f is a map from to 12® ®^ of the 
form f = C + N, where C is a tridiagonal operator satisfying (§, (H and ([^, and where 
the non linear part N is quadratic. Assume that for some m > 6 we have computed an 
approximate zero of f, of the form x = {xq, ..., Xm-i, 0,..., 0,...), and D an approximate 
inverse ofDf^) {x). Consider 


T : 


12 ® ^ 12 ®, 

xe^ X — Af{x), 


where A is defined as in (21). Take M satisfying (32) and (36) and L>0 a computational 
parameter. Then the bound Y defined in Proposition 3.10 satisfies (26) and for all r > 0, 
the bound Z{r) defined in Proposition 3.1!^ satisfies (27). 


Now that we have found bounds Y and Z{r) that satisfy (26) and (27), we must find a 
radius r > 0 such that ||y -I- Z(r)\\s < r in order to apply Theorem 3.1 By definition of 
the norm ll’Hj,, it amounts to find an r > 0 such that, for every A: > 0, the radii polynomial 
Pk{r) satisfies 

Pfc(r) = Yk + Zfe(r) - 4 < 0. 

Note that since we constructed Y and Z in such a way that for every k > M, 


Pm-t-fc — 


m+M ■ 


^m+M 

,s 


a,Ild Z^n-l-k — 


^m+M ■ 


^m+M 

'-’m+k 


it is enough to find an r > 0 such that for all 0 < fc < m -I- M, Pfc(r) < 0. In order to do so, 
we numerically compute, for each 0 < k < m + M, 


Lk = {r > 0 I Pfc(r) < 0}, 
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and 


m+M 

1= n 

k=0 

If I is empty, then the proof fails, and we should try again with some larger parameters 
m and M . If I is non empty, we pick an r € I and check rigorously, using the interval 
arithmetics package INTLAB m, that for all 0 < fc < m + M, Pk{r) < 0, which according 
to Theorem 3.1 proves that T defined in (25) is a contraction on Bs{x,r), thus yielding the 
existence of a unique solution of f{x) = 0 in Bs{x,r). 


4 An example of application 

We present in this Section an example of equation, for which it is possible to apply the 
method developed in this paper. We first explain the link between the equation that we 
study (cf. ( [50| below) and the tridiagonal operator defined in Section]^ Then, we explain 
what are in this example the values of the various constants and parameters of our method. 


Equations of the following form: 

- (2 + cos^)u"{^) + u{^) = 

m'( 0 ) = u' {tt) = 0 , 


(50) 


where 5 is a 27r-periodic even smooth function, fall into the framework developed in Sectionj^ 
Consider indeed the cosine Fourier expansions of u and g\ 

uiO = '^Xk cos(fc^), g{Cs ='^9k cos{kO- 


k^Z 


kGZ 


Then, (50) can be rewritten as f{x) = 0, where 

fo{x) = xo + xi + (J {x * x)q - go, 

and for all /c > 1, 

fk{x) = ^{k- ifxu-i + (1 + ‘2.k‘^)xk + + ^fxk+i +(T[x*x)f^- gk- 

We see that the linear part of (51) is, as in given by 

Pk{x^ — T g^kXk PkXk-\-l, 

with 


def -I ^ def 

Mo = 1, Po = 1, 


and for all /c > 1, 


^k = ^ik-l)^, /ifc = (I + 2fc^) and /3k = ^{k + 1)'^. 


Let us fix some m > 2. With 


Cl =2, C2=3 and d= , 

4 + 4 


(51) 
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we get 


V fc > 1, 




fj-k 


Pk 


1 

/c2 




<C2, 


together with 


y k > m, Cl < 
We now focus on the example when 


i^k 

and 

Afc 


Pk 



gk 


g-k 


< < 5 . 


giO = \ +3cos(0 + ^cos( 2 C), 

so that u{^) = cos(^) is a trivial solution for tr = 0. We are going to use rigorous com¬ 
putations in order to prove the existence of solutions for cr 7 ^ 0 , and to compute these 
solutions. 

Starting from <7 = 0, we first use standard pseudo-arclength continuation techniques 
to numerically get some nontrivial approximate solutions for cr 7 ^ 0. We computed 1250 
different solutions (675 for cr > 0 and 675 for a < 0). See Figure for a diagram summing 
up those computations, where each point represents a solution of (50). 



Figure 2: Branch of solutions of (50). 


Then we use the rigorous computation method described in this paper to prove, for each 
numerical solution, the existence of a true solution in a small neighbourhood of the numerical 
approximation. We keep m = 20 Fourier coefficients for the numerical computation, and 
use M = 20 and the decay rate s = 2 for the proof. The bounds of Lemma 3.12 as well as 
the error on (35) are computed with L = 100. For each numerical solution, the proof is 
successful. The set I defined in Section on which all radii polynomials should be negative 
always contains [4 x 10“^^, 10“^], and we rigorously prove using interval arithmetics that 
they are indeed all negative for r = 10“^°. Hence the assumptions of Theorem 3.1 hold 
and as a consequence, within a ball of radius r = 10 “^° in H'* centered on the numerical 
approximation, there exists a unique solution to (50). Therefore the existence of the solutions 
represented in Figure is rigorously proven, within a margin of error that is too small to 
be depicted. The codes used to perform the proofs can be found in m- 
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Notice that existence of solutions of ( [50| ) could certainly have been obtained in different 
and more classical ways, for example using perturbative methods when cr is close to 0, 
or using a variational approach (that is, considering (50) as the Euler-Lagrange equation 
related to the critical points of a functional), or even using topological tools such as the 
Leray-Schauder theory. The advantage of our method is that it gives us more quantitative 
information than those approaches: indeed it enables to provide more than one solution for 
some values of cr, and, maybe more importantly, it gives a very precise localization of this 
(or these) solution(s) in terms of Fourier coefficients (something that looks very hard to 
obtain with qualitative PDEs methods). 


5 Conclusion and Perspectives 

A first interesting future direction of research would consist in adapting our approach to the 
rigorous computation connecting orbits of ODEs (using spectral methods). For instance, 
we would like to investigate the possibility of combining Hermite spectral methods with our 
approach to compute homoclinic orbits (e.g. see [Tsms]). Since the differential operator in 
frequency space of the Hermite functions is tridiagonal, adapting our method to this class 
of operator could lead to a new rigorous numerical method for connecting orbits. 

It would also be interesting to adapt our method to the case of solutions belonging to 
the sequence space 


i\={x= {Xk)k >0 ■ \\x\\j, = ^ \XkW^ < oo} 

k>0 

for some v > 1. With this choice of Banach space, we could use the fact that il is naturally 
a Banach algebra under discrete convolutions. This could greatly simplify the nonlinear 
analysis. 

Note that assumption (j^ requires the tridiagonal operator to have symmetric ratios be¬ 
tween the diagonal terms and the upper and lower diagonal terms. This is a restriction that 
could hopefully be relaxed. Since many interesting problems involve tridiagonal operators 
with non symmetric ratios (as in the case of differentiation in frequency space of the Hermite 
functions), we believe that this is a promising route to follow. 

Finally, generalizing our approach to problems with block-tridiagonal structures could 
also be a valuable project. 
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